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A problem of determining of the relaxation time distribution in irregular relaxing systems is 
considered from the point of view of the Kotelnikov-Shannon theorem of the telecommunication 
theory. 

The dielectric relaxation in irregular systems like polymers, impure ferroelectrics etc, and magnetic relaxation 
in impure magnetics including spin glasses [lj are usually described with a help of a continuous relaxation time 
distribution function (RTDF). It is desirable to have, along with theoretical models for such distributions, some regular 
technique for determining of them using the measured experimentally dielectric e(w) or magnetic n (oj) spectrum. We 
will discuss for definiteness the former case. The exact formula relating the dielectric losses and RTDF was derived in 
many years ago . However, this formula can not be used directly in practice (see the discussion below). Numerical 
analysis of the proper integral equation was done recently 0|, 0]. Our aim here is to develop an approach giving a 
practical procedure for determining of RTDF. 

We introduce RTDF g(ln(r/ro)) with a help of the relation 



X(M) = / XD(uT)g(ln(T/To))dln(T/T (1) 
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with 

Xd(ut) = = : — , (2) 
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where eo (wt) is the Debye single-relaxator dielectric permittivity, eo and are respectively zero- and high-frequency 
limits of the permittivity, tq is an arbitrary time scale. Basing on rather intuitive grounds, the distribution function 
is often estimated assuming g(ln(T/ro) = Imx(ln(wro). Our aim here is to give a regular procedure for calculation of 
the distribution function (;(ln(r/ro) basing on a measured dielectric spectrum e(uj). 
Let us introduce the logarithmic variables: 

z = \tl(t/t );u)t = VL;t/t = cxp(z),D = cxp(-/i), (3) 
lot = ilexp(z) = exp(z — ^),y = — \xi(ujtq), 

Then equation (1) can be rewritten in the form 

x(y) = + f dx . 9{x) (4) 

Now we discuss some possible approaches to the problem. The integral equation Q can be solved exactly. This 
equation belongs to the class of equations that can be solved using the Fourier transform. Taking the imaginary part 
we have 

X (w = / d\n(T T ) — — . (5) 
Introducing the logarithmic variables we can write 

x " (/i) = p dz g(*)"exp(z) = r dz ,g(z)exp(z- M ) = 1 g{z) 

i-oc l + r» 2 cxp(2z) J_ x l + cxp(2(z- A1 )) 2 7^ cosh(z- M ) W 
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Let us write (JHJ in the standard form 



x "(fx)= / dzK(fi-z)g(z). 

J — CO 

This integral equation can be solved using the Fourier transform: 



(7) 



X y = I dzexp(iyz)x"(z), 



(8) 
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dz exp(iyz)g(z), 



(9) 



K(z - fi) = l/[2cosh(z - fi)]. 
Performing the Fourier transform we obtain an algebraic equation relating the Fourier transformants: 



(10) 



Xy = K v g v . 



Using the known integral 



dx 



cos xy 



cosh ax 2acosh(|^ 



we can write 



K y = 7r/[2cosh(7r?//2)]. 



(11) 



(12) 
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Thus the Fourier transform of the solution reads 



.g y = 27r- 1 X ;'cosh(7ry/2). 



Taking the inverse Fourier transformation 



9(z) 



dy_ 

2tt 



eTq>(~iyz)g y , 



(14) 



we can, in principle, calculate RTDF basing on the measured dielectric losses function x"(a*)- These integrals must 
be computed numerically if the susceptibility is presented as a point numerical function. The trouble with this way 
is the extremely poor convergence of the integral with respect to /i. However, there is another possible way: if we 
can carry out an analytic continuation of the dielectric losses function we can integrate generically. We assume the 
integrals to be uniformly convergent within the domain of integration so that one can change the order of integration 
(an exact analysis must be carried out on the base of the theory of distributions [? ]): 



1 



(z) = — I dy I d/zexp[(/z - z)yi]x" (fi) [exp(7ry/2) + exp(-7n//2)] 



dfix (p) [5{z — 11 + in/ 2) + S(z — fj, — wr/2)] 
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= i[ x "(z + in/2)+ X "(z-m/2)} (15) 

7T 

This solution was derived by Kirkwood Q in a different way. 

Now we discuss possible ways for practical use of this formula. Assuming differentiability of x"( z ) U P to at least 
the second derivative and its slowliness and expanding it in powers of «7r/2, we can give the approximate form of the 
solution: 



7(*) « - 

7T 



^ d 2 X "(z) 
8 dz 2 
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This approximate formula is asymptotically exact for very slowly- varying dielectric functions. Formula 1151 
can be basically used in a general case, but ambiguity of the analytic continuation makes it problematic. There are 
two main reason for it. Firstly, results of the analytic continuation essentially depend on analytic properties of the 
function, but they cannot be established exactly basing on a function determined on the finite segment of the real axis. 
Secondly, measuring of the dielectric losses function yields really a discrete set of points that presents an additional 
ambiguity factor. This means that analytic properties must be determined independently. A similar problem exists in 
the communication theory. Discretization of a signal and impossibility to transfer a signal with unlimited spectrum 
through the channel make it necessary to apply the approach developed in the communication theory by Kotelnikov 
and Shannon ;|| (basically similar approach was really known earlier in the mathematical theory of approximation). 
Following to the spirit of this approach, we assume that x" ( z ) is an entire function of finite degree Q : 



lnmaxi z |_ r |y"(z)I , , 

lim N , V n = a<+oo (17) 

But according to the Paley- Wiener theorem |(|, such a function has a bounded spectrum. What does it mean from 
the physical point of view in our case? The spectrum of x"i z ) characterizes the time evolution of the polarization; 
therefore, boundness of the spectrum means that a reaction of the polarization on the field pulse decays for some 
finite time. Surely, this decay time is theoretically infinite in glasses, but it is essentially restricted by noise and other 
factors in a real experiment. That is why we have right to assume x"( z ) to be an entire function of finite degree. 

Let x"{ z ) be a function with the spectrum g y given by El lying in the band —j3 < y < (3. Then the following 
interpolating formula holds for x"( z ) at a > /3: 



... , 7t sin a (z — k — ) , . 

This series converges uniformly to the entire analytic function and admits an analytic continuation. Substituting 111 5(1 
into iJTSJ we come to the final formula: 
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sin a (z + i% — k~\ sin a (z — — k—) 
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(19) 



Thus, we have shown that for the case of the finite-time response functions, the "naive" solution g(z) oc x"( z ) is 
exact within the class of the entire analytic functions. Notwithstanding the fact that a real measurement gives only 
finite-time response because of the noise effect, it would be interesting to generalize this result to the case functions, 
which have an infinite spectrum with power-law decreasing. We hope to return to this problem in the future. 
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